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a^ : ABSTRACT 

The possibility of detecting supersymmetric dark matter is examined 

^ . within the framework of the minimal supergravity model (MSGM), where the 

Zi is the LSP for almost the entire parameter space. A brief discussion is 



o 



given of experimental strategies for detecting dark matter. The relic density is 
constrained to obey 0.10 < O^ h <0.35, consistent with COBE data. Ex- 

o 

pected event rates for an array of possible terrestial detectors ( He, CaF2, Ge, 

GaAs, Nal and Pb) are examined. In general, detectors relying on coherrent 

f^ , ^1 -nucleus scattering are more sensitive than detectors relying on incoherrent 
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> 



in 



(spin-dependent) scattering. The dependence of the event rates as a function 

of the SUSY parameters are described. The detectors are generally most sen- 

, sitive to the small UIq and small ITlq and large tan/3 part of the parameter 

Qs . space. The current b — ^ S -\- 'J decay rate eliminates regions of large event 

(-H , rates for |U > 0, but allows large event rates to still occur for // < 0. MSGM 

CIh! models that also possess SU(5)-type proton decay generally predict event rates 

O ' below the expected sensitivity of current dark matter detectors. 

0-) ■ 

1. Introduction 



^ 



/\ ' If the SUSY models currently being examined are correct, then supersymmetry 

c^ I will be discovered at the LHC in the year 2005 or possibly even at an upgraded 

version of the Tevatron (e.g. the DiTevatron) in the year 2000. However, high 
energy colliders may not shed further light until then. Thus it is of interest to 
look at other phenomena which supersymmetry might effect, e.g. dark matter, 
proton decay, the b ^ s + •y decay etc. Each of these restricts the parameter 
space of supersymmetry, and so by combining the constraints, one can get sharper 
predictions of what to expect at colliders. We will discuss here the question of 
detecting SUSY dark matter, and how bounds on other processes affect dark matter 
searches. There is a warning however one should make concerning such analyses: 
In applying SUSY to dark matter searches, one is making additional assumptions 
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(e.g., cosmo logical assumptions) not made in dealing with accelerator phenomena. 
We will test the sensitivity of final results to some of these extra assumptions, but 
some caution is needed in interpreting the theoretical predictions. 

2. Dark Matter 

There is much astronomical evidence that more than 90% of the total mass of 
the Galaxy and perhaps of the universe is made up of dark matter of an unknown 
type. In the Galaxy this can be seen from rotation curves of luminous matter (Fig. 
1). The circular velocity Vdr does not fall with r beyond the optical radius. Similar 



Fig. 1. Estimated mean spherical density of dark matter in the Galaxy^. 

effects are seen in other galaxies (Fig. 2). In the vicinity of the sun, the mean 
density of dark matter (DM) is estimated as 

Pdm = 0.3 GeV/cm^ (1) 

(about 10'^ 'J universe) and assuming a Maxwell velocity distribution, the DM has a 
velocity relative to the solar system of 

VDM = 320A;m/s (2) 



i.e. vdm/c ^ 10" 



Fig. 2 Vcir for a number of galaxies showing Vdr remains approximately constant 
well beyond the optical radius^. 

There are a large number of candidates for dark matter, both from astronomy 
and particle physics. Supersymmetric models with R-parity offer two candidates: 
the lightest neutralino Zi and the sneutrino, z/. However, in supergravity models 
almost always the Zi is the lightest supersymmetric particle (LSP) and hence is 
absolutely stable. Thus the relic Zi, left over from the big bang would be the dark 
matter. The dynamics is then well fixed and we will deal with this case exclusively 
here. 

COBE data suggests that DM is a mix of cold dark matter (CDM) [which 
we are assuming here to be the Zi] and hot dark matter (HDM) [possibly massive 
neutrinos] in the ratio of about 2:1. There can also be baryonic dark matter (B) 
but nucleosynthesis analyses limit this to ~ 10%. Thus if we define Vti = Pi/pc, 
where pi is the mass density of type i and pc = 3H'^ / {SttG n) is the critical mass 
density to close the universe (H=Hubble constant and GAr=Newtonian gravitational 
constant), the infiationary scenario requires EO^ = 1 and hence a reasonable mix 
of matter is 

n^, - 0-6; flHDM ^ 0.3; Qb ^ 0.1 (3) 

What can be calculated theoretically is O^ /i^ where h=H/(100 km/s Mpc). 
Astronomical observations give a range of values for h i.e. h = 0.5 — 0.75. Hence 



O^^/i^ ^0.10-0.35 (4) 

This bound strongly restricts the SUSY parameter space. We will discuss below 
how sensitive our results are to the endpoints of this bound on Q^; ^^• 

3. Detection Strategies 

The solar system is presumably being bombarded with Zi particles moving 
with velocity < v^ >= 320 km/s. Two strategies have been proposed for their 
detection. 

3.1 Indirect Detection 

The Zi impinging on the sun can be relatively easily captured since the escape 
velocity at the surface of the sun is 618 km/s= z/^ . Once captured the Zi will be 
slowed down inside the sun and gradually fall to the center. There they accumulate 
and can annihilate giving rise to neutrinos via e.g. Zi + Zi —^b + b,....—^u^ + X. 
Since m^ =(10-100) GeV, high energy neutrinos coming from the sun from this 
process would be a striking signal that could be observed on Earth by a neutrino 
telescope. Calculations indicate that one would need a telescope of area > Ikm'^ 
to cover the SUSY parameter space^, and telescopes of this size are currently being 
built. 

3.2 Direct Detection 

A direct approach to see the incoming Zi is to detect their scattering by quarks 
in nuclei of a terrestrial detector: Zi + q ^' Zi + q. Two types of detectors being 
considered are the following. 

3.21 Low Temperature Detectors 

Since m^ — (10-100) GeV, the recoil energy to a nucleus that has been struck is 
AE ^ {v^ /c)^(to^ c^) — (10-100) KeV. This is of a size to produce phonons (heat) 
or ionization in the lattice of the detector. The temperature rise AT is AT=AE/VC 
where C is the specific heat and V the volume. Since at low temperature C~T"^ 
one can enhance AT by reducing the temperature but one is also limited by V not 



being too large. An optimum set of parameters is to have T in the mK range and 
a detector of mass of ~lkg. Detectors of this type are currently being built. 

3.22 Superconducting Detectors 

Here superconducting granules are suspended in a dielectric carrier in the pres- 
ence of a magnetic field (Fig. 3.) The superconductors are put into a metastable 



Fig. 3 Schematic diagram of a superconducting detector^. 

state and hence the magnetic field is excluded from the granules by the Meissner 
effect. When a Zi strikes a granule, the deposited energy triggers the transition to 
the normal state, and the magnetic flux movement then produces a signal in the 
pickup coil. The characteristic size of such detectors are also about 1 kg. 

Background for these DM detectors include cosmic ray niuons and natural 
radioactivity. The present sensitivity expected is 0.1 events/kg da, and this might 
be improved at a later date (i.e. by going underground) to 0.01 events/kg da. 

4. Dynamical Model 

To calculate the event rates expected at terrestial detectors, we need to calcu- 
late two items: (1) the relic density of the Zi, in order to make sure that O^ h'^ 
lies in the range of Eq. (4), and (2) the Zi-nucleus cross section to obtain the 
expected event rate for a given detector. The bounds on O^ h'^ signiflcantly limits 
the SUSY parameter space and hence strongly affects the event rates obtained from 
the Zi-nucleus cross sections. 

In order to carry out the above calculations one needs a dynamical model and 
we use here models based on supergravity grand unification"^. These models have 



the advantages of (i) being consistent with the LEP data on unification of couphngs 
at Mq — 10^^ GeV, (ii) generating spontaneous breaking of supersymmetry at 
Mgr (in the "hidden" sector), (iii) generating naturaUy spontaneous breaking of 
SU(2)xU(l) at the electroweak scale by radiative corrections, and (iv) having all 
new SUSY phenomena described by only 4 new extra parameters and one sign. 

The minimal supergravity model (MSGM) is characterized at the Gut scale by 
a superpotential 

Mg 

where Wy is the cubic Yukawa couplings and W*-^-* is the quartic non-renormalizable 
couplings (possibly leading to proton decay). In addition there is a soft supersym- 
metry breaking effective potential 

VsB = vr?^Za^Za + {A^Wy + Bo^oHiH2 + h.c.) (6) 

a 

where {za} are the scalar fields, and a universal gaugino mass term L^^^^ — 
— rMi/2A"A" where A" are the gaugino fields. The scalar mass mo (and cubic soft 
breaking constant Aq) are universal provided the agent of supersymmetry breaking 
in the hidden sector (e.g. the super Higgs field) communicates with the physical 
sector in a fiavor independent way. This is automatically the case for contribu- 
tions arising from the effective potential (since there the only communication is 
gravitational) and will be true in general if the couplings of the superHiggs to the 
physical fields in the Kahler potential is also fiavor independent. Eq. (6) guarantees 
a natural suppresion of unwanted FCNI. 

Using the renormalization group equations (RGE) one obtains at the elec- 
troweak scale the conditions for SU(2)xU(l) breaking 

-Ml = -^^ + ^^ , ^^ -; sin2(3 = . ^ ^ 7 

where tan [3 =< H2 > / < Hi > and m'jj, are the Higgs running masses including 
1-loop corrections. One can then eliminate Bq and /Iq leaving four parameters toq, 
mi/2. At (the t-quark A parameter at the electroweak scale), tan/? and the sign of 
Ho to determine all the masses and interactions of the 32 new SUSY particles. One 
therefore expects a large number of relations holding between the SUSY masses. If 
one limits the parameter space so that (i) experimental mass bounds of LEP and the 



Tevatron are obeyed, (ii) mo, rrig < 1 TeV, rrig the gluino mass (so that no extreme 
fine tuning of parameters occurs), and (iii) radiative breaking of SU(2)xU(l) occurs 
at the electroweak scale, then the foUowing "scahng" relations hold throughout most 
of the parameter space^ : 2m^^ = m^^ = m^y^ ^ (^"i)"^;?; "^z ~ "^^4 ~ "^VK2 ^^ 
7712; '■> ^^^ niHO = niAo = mH± » TTih- Here the Wi are the charginos, the Zi are 
the neutralinos, h^ and H^ are the CP even Higgs, A'^ is the CP odd Higgs and H^ 
is the charged Higgs. In addition one always has tan/? > 1. 

5. Calculation of Zi Relic Density 

As discussed above, Eq. (4) puts a significant constraint on the allowed region 
of parameter space. Since this band is relatively narrow, it is important to include 
major effects in calculating the relic density. 

R parity makes the Zi produced in the early universe absolutely stable. How- 



Fig. 4 Annihilation diagrams of the Zi in the early universe. 

ever, they can annihilate in pairs and the main annihilation diagrams are shown 
in Fig. 4. The calculation of the mass density of the Zi remaining at the present 
time proceeds as follows^: Initially the Zi are in thermal equilibrium with the 
background and the reactions of Fig. 4 go forward and backward. However, when 
the annihilation rate falls below the expansion rate of the universe, "freezeout" 
occurs at temperature Tj, and the Zi disconnect from the background. The Zi 
then continue to annihilate and the amount left at present time is^ 



n.h' = 2.4 X 10— ( ^ {(^Xtt^ (8) 



-'— -"(^) (l#) ^, 



where Nf is the effective number of degrees of freedom. 



J{xf) = < av > dx; x = kT/m^^, (9) 

Jo 

and < (jv > is the thermal average (cr = annihilation cross section, v — relative 
velocity): 

/'OO /'OO 

< av >= / dvv^ (av) Exp{—v^ / Ax) / dvv^ Exp{—v'^ / Ax) (10) 

Jo Jo 

In general, the annihilation occurs non-relativistically, i.e. x/ ~ ^. However, this 
does not mean one can always make a non-relativistic expansion, av = a + bv^ + ..., 
in performing the calculation of J(xj). As has been pointed out^, such an expansion 
fails near an s-channel pole. For example for the h pole one has 

{s — Am'2 )/w| 
av — A- TT-^ ^-V (11) 

where r^=0(MeV) is the h boson width and A is a constant. Since Fh is small one 
must treat the pole more carefully. The danger of not doing so is shown in Fig. 5^. 



Fig. 5. ^approx is the Zi relic density calculated in the av = a + bv'^ approximation 
and O is the relic density calculated rigorously using Eq. (11). The h and Z poles 
occur at the points where the curve descends through zero. 

Note the long tail where an error of a factor of ~ 2 can be made well past the h and 
Z poles. One finds, for example, that ^approx has an error of > 25% for 65% of the 
parameter points where rUg < 450 GeV, while ^approx is a good approximation for 
~ 100% of the points for nig > 450 GeV. The reason for this is the scaling relations 
discussed at the end of Sec. 4. Since 2m ^ > (| — |)?Tig, when rrig > 450 GeV, the 
Zi has passed both the h and Z poles (where 2m^ ~ rrih or Mz). (Recall that in 



the MSGM one has rrih < 130 GeV.) However, for hghter nig the pole effects are 
very important, since one is almost always near an h or a Z pole. 

6. Detector Rates 

The dark matter detectors discussed in Sec. 3.2 detect the Zi from their 
scattering by quarks in the nucleus. The basic diagrams are given in Fig. 6. They 



Fig. 6 Zi-quark scattering diagrams, q is the squark. 

are mainly crossed diagrams to the relic annihilation diagrams of Fig. 4. Thus to 
a rough approximation, when the relic density of the Zi is small (i.e. there is a 
large annihilation cross section) the number of scattering events will be large. This 
makes the results somewhat sensitive to the lower bound, O^ /i^ > 0.1, that we 
have chosen, and we will discuss this sensitivity below in Sec. 8. 

Calculations show that it is possible to represent the Zi-q scattering by an 
effective Lagrangian^: 

Leff = {xiYl''Xi)qi^{A,PL + BqPR)q + (xiXi)C,m,Qg (12) 

where Xi(x) is the Zi field, q(x) the quark field, and ruq its mass and Pr^l = 
^(1 ±7^). The coefficients Aq, Bq arise from the Z t-channel and q s-channel poles 
while Cq arises from the h and H t-channel and q s-channel poles. 

The first term of Eq. (12) leads to spin dependent incoherrent scattering by the 
quarks in the nucleus, while the second term leads to coherrent scattering, where 
the masses of all the quarks (and hence nucleons) approximately add coherrently. 
Thus the coherrent amplitude will contain a factor M^ larger then the incoherrent 
amplitude, where M^ is the nucleus mass. 

In general, the Zi field is a mix of higgsinos and gauginos: 



Xi = niWs + n2B + n3H2 + n4:Hi (13) 

where W3 and B are the W^ and U(l) [Bino] gauginos, and Hi 2 are Higgsinos. The 
Hi can be computed in terms of the basic SUSY parameters niQ, mi, At and tan/3. 
Over most of the aUowed part of the parameter space one finds 

n2>ni,n4»ns (14) 

In SUSY theory, nih is smaU i.e. nih < (120 — 130) GeV, and niH^ » ?ti^2. In 
spite of this it is remarkable that the H contribution is important in the coherrent 
amphtude^^. The reason for this is the following. One finds for C„ the result 



c„- »' 



' 4Afw2 



cosa P'h I I sing ^H 

sinf3 ^h K I ) sinp rn^ 

sing Fh ( | cosot Fh 

cosf3 mT I I ^°^f^ m^r 



u — quark , , 

d — quark 



where a is the rotation angle needed to diagonalize the 2x2 h-H mass matrix, 
Fh = {ni—n2tan6w)ij^2,cosa~n4^sina) and F// = (?ii— ?i2tan6'v7)(n3sma— ?i4Cosa). 
In calculating the h-H mass matrix, one must include as is well-know, the one loop 
corrections due to the fact that rrit is large. One finds then, for the allowed part 
of the parameter space that a is generally small i.e. a ~ ^. This result combined 
with Eq. (14) shows that generally cosaFn » sinaFh, which can overcome the 
reduction in the d-quark amplitude due to the largeness of mn- One finds, in 
general, as one varies over the parameter space that the H contribution to Cq can 
vary from ^ to 10 times the h contribution for d-quarks, but is generally a small 
correction for u-quarks. 

7. Detector Event Rates 

The total event rate expected for a given dark matter detector can be written 
in the following form^ 

R-iR...^R...]i.r..m:^nJ:i^f4^ (16) 



where 



0.3GeVcm ■^ 320/cm/s kg da 



4m^M%M% 



p — ^1 " ^ I /I 1-=! n 7^1 
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R^nc = JTpj^^^X^J{J+l)\A,nc\' (17a) 

[Mat + 7712^-^ 

Here Acoh ~ Cq, Ai^c ^ Bq — Aq, J is the nuclear spin and A < A^| J|A^ > 
= < A^|E5'i|A^ > where the sum is over the spins of all nucleons in nucleus N. 
Note that for large Mjv, Rcoh increases as M^ while Rinc decreases as l/M^. This 
additional M^ factor in Rcoh is as expected from the discussion following Eq. (12). 

We have examined the expected detector event rates for detectors made from 
the following nuclei^^: ^He, "^^Ca ^^^2, ^^Ge + ^^Ge (equal mix of isotopes), 
'''^Ga ^^As, ^^Na ^^^I, and ^'^^Pb. Both ^^F and ^He have strong spin interactions, 
while Ge, I and Pb are increasingly heavy and hence have increasingly strong co- 
herrent scattering. 

The parameter space studied was 

lOOGeV < mo, nig < ITeV; -6 < At/nio < 6; 2 < tanP < 20 (18) 

and the mesh used was Amo=100 GeV, Amg=25 GeV, AAt/niQ = 0.5, A{tan(3)=2, 
4, and the t quark mass was set at mt=l67 GeV. 

The dependence of the event rates on the SUSY parameters is generally quite 
complicated. However, it is possible to understand these dependences in a qualita- 
tive way. Fig. 7 shows that event rates decrease rapidly with rrig. This follows from 
the fact that the Zi becomes increasingly Bino as rrig increases i.e. n2 of Eq. (13) 
grows and hence ns and n^ shrinks, making the interference between the gaugino 
and Higgsino in F^ and Fh (needed for coherrent scattering) become small. CaF2 
has the strongest spin dependent forces while Pb is the heaviest of the detectors 
chosen. One sees from this that Rcoh is significantly larger than Rinc^ a general 
feature. 

Fig. 7 also shows that R increases with tan/?, a feature which can be seen in 
detail in Fig. 8. This behavior follows from the l/cos/3 ~tan/3 factor for d-quark 
scattering in Eq. (15). One sees again the scaling of R with Mn i.e. the ^^Na^^^I 
curves lie above the ^"^Ge+^^Ge curves. The dependence of R on ttiq is somewhat 
complicated. One expects R to decrease with toq since m? increases with m^ and 
hence the squark s-channel contribution of Fig. 6. shrinks. In addition, however, 
the iJ? determined by the radiative breaking condition Eq. (7) increases as mg 
increases, making the Zi increasingly more Bino like which further reduces Rcoh- 



11 



This decrease in R with increasing mo shown in Fig. 9. 



Fig. 7 Event rate as a function of nig for Pb (sohd) and Ca F2 (dashed) detectors. 
The upperhne in each pair is for tan/3=20 and lower hne for tan/?=6. The curves 
are for At/mo=1.5, mo=100 GeV, fx > 0. 



Fig. 8 Event rates for Nal and Ge detectors vs. tan/? for mg=275 GeV. The dot dash 
curve is for At/mo=1.0, ?bo=200 GeV, the dashed curve for At/mo=0.5, ?bo=300 
GeV, and the sohd curve for At/mo=0.0, rMo=200 Gev, fx > 0. The upper curve of 
each pair is for Nal, the lower for Ge. The At, nio parameters were chosen so that 
O^ /i^ is approximately the same in each case of a fixed tan/?. 
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The Zx mass is an increasing function of nig and as to 5 increases, the rehc density 



Fig. 9. Event rate as a function of toq for mg=300 GeV, At/rMo=0.5, tan/3=8, 
H > 0. The soUd curves from bottom to top are for Ge, Nal and Pb and the dashed 
curve is for CaF2. 

increases (i.e. the annihilation cross section in the early universe decreases). The 
upper bound on O^ /i^ then leads to an upper bound on nig. In general this upper 
bound on rrig of about 750 GeV. This effect is shown in Fig. 10. 



Fig. 10 {mg)max vs. tan/?, (U > 0. The curves, from bottom to top are for 
At/niQ = -0.5, 0.5, 1.0, 2.0, 3.0 and 3.5. 

One may scan the entire parameter space to obtain the maximum and minimum 
event rates as a function of Af. These are shown in Fig. 11 for yu > and Fig. 12 
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for |U < for the domain 2 < tan/3 < 20. As expected from Fig. 8, the maximum 
event rates occur for tan/?=20. However, the minimum rates occur at different tan/? 
for different Af. Current dark matter detectors can achieve a sensitivity of R~0.1 



Fig. 11 Maximum and minimum event rates as a function of At for // < for Pb 
(sohd) and CaF2 (dashed) detectors. 



Fig. 12 Same as Fig. 11 for n > 0. 

event/kg da with perhaps a factor of 10 improvement in future sensitivity. One sees 
that only part of the parameter space i.e. the region with relativity large tan/3 will 
be accessible. The detectors with large Rcoh (e.g. Pb) are generally considerably 
more sensitive than those with large Rinc (e.g. CaF2). 
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8. Sensitivity To Bounds on O^ /i^ 

In the preceeding discussion, we have assumed the bounds of Eq. (4) for O^ K^ . 
The results have some sensitivity to the choice of endpoints and we examine this 
here. 

As discussed in Sec. 6, small Vt^ h"^ generally leads to large event rates, and 
R generally rises rapidly as O^ h? decreases near its lower bound. The largest R, 
however, occurs for small rrig (see e.g. Fig. 7), and hence by the scaling relations 
discussed at the end of Sec. 4, for small rriy^ . However, there also exist cuts, e.g. 
rriy^ > 45 GeV, required by LEP data which forbids niyy from getting too small. 
Thus one has a sharply rising function R hitting an experimental constraint on the 
parameter space. This is the origin of the sharp peaks in the maximum event rate 
curves of Figs. 11 and 12, i.e. the maximum event rate gets quite large or not 
depending on whether or not the parameter point passes the experimental cut. 

The sensitivity of this effect is seen in Fig. 13 (for Pb) and Fig. 14 (for CaF2) 
for the maximum event rates when one increases the minimum value of O ^ /i^ 



Fig. 13 Maximum event rate as a function of At /mo for ^ <0 for Pb detector. Solid 
curve is for Q.^ h? >0.10 and dot-dash for Q.^ h? >0.15. 

Zi Zi 

from 0.10 to 0.15. The sharp peaks get clipped off when the lower bound on (O^ h?) 
is increased. There still remains, however, a sizable portion of the parameter space 
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where R is large. 



Fig. 14 Maximum event rate as a function of At/nio for n <0 for CaF2 detector. 
Dashed curve is for O^ /i^ >0.10 and sohd curve for O^ /i^ >0.15. 

Zi Zi 

The minimum event rates are sensitive to the upper bound chosen for O^ /i^. 
This can be seen in Fig. 15. The minimum event rates increases by more than 
a factor of 10 as (O^ h'^)max is reduced from 0.35 to 0.20. One notes that the 
inflationary scenario favors a small Hubble constant, i.e. h=0.5, so that the pre- 



Fig. 15 Minimum event rates as a function of (O^ h?)max for Pb detector (solid 
curve) and CaF2 detector (dashed curve), // <0. 

dieted age of the universe not be inconsistent with estimated ages of globular star 
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clusters. Even if one assumed O^ = 1 (i.e. no hot dark matter and negligible 
baryonic dark matter) this would imply O^ /i^=0.25. Thus the larger minimum 
event rates of Fig. 15 may possibly be the correct choice. 

The reason the minimum event rate increases with decreasing (O^ h'^)max is 
that the minimum rates occur for the maximum values of nig. (As seen in Fig. 7 
the event rate drops with increasing rrig.) Further, as discussed in Sec. 7, Fig. 10, 
rrig possesses a maximum value because of the upper bound on O^ h'^. Fig. 16 
shows the dependence of of {rrig) max on the maximum value of O^ h'^. As 



Fig. 16 Maximum value of rrig as a function of the upper bound on O^ /i^ for n <0, 
At/mo=0.5, tan/?=6. Results are insensitive to the values of At and tan/?. 

(0/i^)m,ax is reduced from 0.35 to 0.20, {rng)max drops from 750 GeV to 400 GeV. 
The inflationary scenario thus favors smaller values of rrig. Such low rrig implies 
that the gluino could be detected at suggested energy upgrades of the Tevatron 
(e.g. at the DiTevatron^^). 

9. The b^ s + 7 Decay 

Recently, the CLEO Collaboration have measured the branching ratio for the 
inclusive decay B— > Xg + 7: 

BR {B ^ Xs+-f) = (2.32 ± 0.51 ± 0.29 ± 0.32) x 10"^ (19) 

where the first error is statistical and the last two errors systematic. Combining 
all errors in Gaussian quadrature, one has in the spectator approximation that 
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BR(b-^s + 7) = (2.32±0.67)xl0-^. The b-^s + 7 decay is of particular interest 
in that it begins at the loop level as it is a FCNC process. This means that SM 
and new physics effects enter at the same loop level and one could expect large 
[i.e. 0(1)] new physics corrections to the SM predictions. Thus the b^s + 7 decay 
is an excellent process for detecting new physics. We will investigate what effects 
the current experimental value has on the SUSY parameter space and hence on the 
expected dark matter detection rates. 

The elementary diagrams at the electroweak scale n = My^ are shown in Fig. 
17. The W~ — t intermediate state is the Standard Model contribution, while the 



Fig. 17 Elementary penguin diagrams for b-^s + 7 decay at the electroweak scale 
/J. = M\Y. Only the third generation quarks and squarks make a significant contri- 
bution. 

H~ — t and W~ —t represent additional SUSY contributions. The interactions can be 
represented by an effective Lagrangian for a transition magnetic dipole interaction^'^: 

Le// = GFVti,VZmi,A^SLcr^''bRF^,, (20) 

where the coefficient A^ can be evaluated in terms of the basic parameters too, nig, 
At, tan/3. In order to calculate the decay rate, however, one must use the RG equa- 
tions to evaluate the amplitude at the b scale fx ~ nib. This causes operator mixing 
with the gluonic transition magnetic moment operator SLO'^'^T°'hRG°ij.j^ (where G?^,^, 
a=1...8 is the gluon field strength) and six four quark operators. It is convenient to 
consider the ratio 

^_ BR{h^s + ^) ^ BR{B^Xs + ^) 

BR{b -^ c + e + Ue) BR{B -^ X^ + e + p^) ^ ' 
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since poorly known CKM matrix elements and nib factors cancel out in the ratio. 
One can recover the 6 ^ s + 7 rate then from the experimental number of the charm 
semi-leptonic rate: BR(B^ Xc + e + z/e)=(10.7±0.5)%. To leading order (LO) QCD 
the value of R is^"^'^^ 



^_6a\VZVtt\' 



r]^Cj{Mw) - Uv^ -V^)C8{Mw)+C2{Mw) 



TT iKfcP 



Hz) l-^^'^^fiz) 



(22) 



where 77 — 03 (M^) /eta (rrib) =0.548, z=mc/mb=0.316±0.013, I(z) is a phase space 
factor for the b— ^cez/e decay, Cj{C8) are the Wilson coefficients for the photonic 
(gluonic) magnetic penguin operators and C2 comes from the operator mixing with 
the 4-quark operators. 

There are a number of theoretical uncertainties in the above calculation which 
can be sumarized as follows: (i) Errors in input parameters i.e. as{Mz), mb/mc, 
CKM factor, BR(B -^ XcCUe)', (ii) Errors in the spectator approximation; (iii) There 
are large NLO (next to leading order) QCD corrections. This can be estimated by 
letting n vary from mb/2 to 2 nib and are seen to be about ±25%; (iv) Heavy 
mass threshold corrections in running the RGE from the t quark/squark, H~, W 
threshold^^'^^ down to nib These are about 15% for the t-quark and estimated to 
be ±15% for the SUSY thresholds'^. Thus current theory has an overall error of 
about ±30%. 

The current CLEO measurement of the b— *> s + 7 rate has a significant effect 
on the expected dark matter detector counting rates. Fig. 18 shows the expected 
b^ s ± 7 rate for a characterisitic choice of parameters'^. The BR(b-^ ■§ + 7) is 
increased when fj, and At have the same sign relative to the value when n and At 
have the opposite sign. (This effect comes from the W — t diagram of Fig. 17). 
Thus regions where n and At have the same sign can exceed the CLEO measured 
value and such regions of parameter space are then experimentally eliminated. One 
sees from Fig. 18 also that the BR(b-^ ■§ ±7) is largest when mo and nig are small, 
which we also saw is the region when the dark matter counting rate R is largest. 
Thus one expects that the maximum values of R get eliminated when /j. and At 
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have the same sign. This may be seen in Fig. 19 (for // < 0) and Fig. 20 (for 



Fig. 18 BR(b^ 'S + 7) as a function of m^^^ for tan/?=5.0, \At/mo\=0.5, mt=165 
GeV, a^^ =24.11. Graphs (a) and (b) are for At < 0,(c) and (d) for At > 0, while 
(a) and (c) are for (U > and (b) and (d) for n < 0. 

H > 0). There we have plotted Rmux for the Pb detector without the b-^ s + 7 
condition and Rmox with parameter points excluded when the predicted b^ s + 7 



Fig. 19 Rmux vs At/niQ for Pb detectors for fx <0. The solid line is the expected 
rate without b-^ s + 7 constraint and the dashed line is the rate with parameter 
points excluded where the predicted b-^ s + 7 rate lies outside 95% C.L. bound of 
the experimental value of Eq. (19). 
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rate exceeds the 95% C.L. of Eq. (19). We see that when fx and At have the same 
sign, the maximum event rate drops sharply, weU below what could be observable 



Fig. 20 Same as Fig. 19 for n >0. 

in the forseeable future. However, the COBE bounds of Eq. (4) turn out to allow 
mostly At >0. Thus the major effect of the CLEO measurement of b^ s + 7 is for 
fj. >0 where for most of the parameter space, the event rate will be very small and 
hence unobservable. However, for fi <0, the b^ s + 7 measurement does not effect 
the expected rates very much, and large dark matter event rates are still possible. 

10. Proton Decay 

The preceeding discussion has been for a generic supergravity Gut model de- 
scribed in Sec. 3 for the parameter domain of Eq. (18). Results are generally 
independent of the Gut group and Gut physics provided that Gut threshold effects 
are not so strong that they prevent grand unification from occuring at Mq ~ 10^^ 
GeV. 

Proton decay is characteristic of all supergravity Gut models except for the 
flipped SU(5) model^^. Further supersymmetry generally implies a unique dominant 
decay mode: 

p^u + K+ (23) 

Thus the observation of this decay would not only indicate the validity of grand 
uniflcation but also of supersymmetry. One can suppress this decay by specially 
tailoring the form of the Gut Higgs sector, but this generally requires some awkward 
fine tuning. We consider here "SU(5)-type" proton decay^^ which arises via the 
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exchange of a superheavy color triplet Higgsino Hs, Fig. 21. (This can happen in 
SU(5), 0(10), Eq etc. Gut groups.) The decay rate can be written as, 

r{p^u + K+) = '^\B\^ (24) 

where B is the loop amplitude. The current experimental bound^^. 



Fig. 21 Example of proton decay diagram for p ^ u + K~^ for supersymmetric 
grand unification. The H^ vertices violate baryon and lepton number. 

T {p ^ u + K) >1.0xlO^^ y, then leads to a bound on B 

B ~ 100 (MhJMg) GeV-^; Mq = 2 x lO^^GeV (25) 

We restrict Mh^ to obey Mhq/Mq < 10 in the following so that the Gut scale 
be disjoint from the Planck scale. (For larger Mho, one might expect large Planck 
physics corrections to enter, the nature of which are not known.) B can be expressed 
in terms of m^^, rriq etc. and hence by the RG equations in terms of the four basic 
parameters mo, m^. At, tan/? and the sign of //. Thus the condition on B is a 
constraint on the parameter space. 

The second generation dominates the loop of Fig. 21 and to a rough approxi- 
mation on has 

^"' ""^ X lO^Gey-i (26) 



azsin2l3 m^ 

where m? = ttiq + 0.6m~. Thus the proton decay bounds of Eq. (25) imply small 
rrig, large tuq e.g. mo > rrig, and small tan/3 i.e. tan/3 <10. One may satisfy both 
the COBE constraints on Vt^ K^ of Eq. (4) simultaneously with the proton decay 
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constraints of Eq. (25) even though mo must be large (which usually is a region 
of small relic Zx annihilation). This is possible since, as discussed in Sec. 5, when 
rrig <450 GeV, (the region also required by Eq. (25) for proton decay) large relic Zx 
annihilation can occur due to the presence of h or Z poles. In fact, one finds that in 
the vacinity of these poles, when too is small (e.g. toq < wig) too much annihilation 
occurs [i.e. O^ /i^=O(10~^)] and one must increase ttiq to satisfy the lower bound of 
Eq. (4). These are then the domains that also satisfy the proton decay constraint. 
Thus in order to find the parameter space region which simultaneously satisfies the 
dark matter and proton decay constraints, it is essential to treat the calculation of 
Zx relic density in the accurate fashion discussed in Sec. 5. One finds then that the 
parameter points satisfying the simultaneous constraints require 

m-g < 375GeV; mo > 500GeV; tan(3 < 10 

0.0 < At/mo < 0.5 (27) 

for a t quark of mass 165 GeV. 

One may next ask whether the CLEO measurement of the b— *> s + 7 decay 
effects this result. One finds, however, that 95% of the parameter points which 
simultaneously satisfy the dark matter and proton decay constraints, predict a 
b— * s + 7 branching ratio in the LO that is within the 90% C.L. of the experimental 
value of Eq. (19). Thus, at the current experimental and theoretical accuracy, the 
b^ s + 7 decay does not effect the proton decay predictions. 

However, the proton decay constraint does effect the expected dark matter 
event rates, and the maximum event rates are significantly reduced since tan/? is 
small and mo is large. Parameter points which simultaneously satisfy COBE, the 
b^ s + 7 and proton decay constraints lead to event rates of size R=O(10~"^- 
10~^) events/kg da. Thus if the next round of proton decay experiments (Super 
Kamiokande, ICARUS) were to actually detect proton decay, the present theory 
implies that Zx dark matter is beyond the ability of current detectors to discover 
by direct detection. 

11. Conclusions 

Studying non high energy accelerator phenomena such as dark matter, proton 
decay, the b— ;► s + 7 decay is useful in limiting the SUSY parameter space. We have 
examined within the framework of the minimal supergravity model (MSGM) the 
ability of direct detection of dark matter when the relic density obeys Eq. (4). 
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(i) The detectors which are most sensitive to coherrent Zx scattering (e.g. the 
Pb detector) are better than the detectors most sensitive to incoherrent (spin- 
dependent) scattering, and the heavier the nucleus, the more sensitive the 
detector. 

(ii) With a future sensitivity of R > 0.01 events/kg da, a reasonable amount, though 
not all, of the parameter space will be accessible to dark matter detectors. 

(iii) Raising the lower bound on O^ /i^ decreases the maximum event rate (Fig. 
13, 14), and lowering the upper bound increases the minimum event rate (Fig. 
15). The upper bound, O^ /i^=0.35, also determines an upper bound on rag of 
750 GeV. Thus if this upper bound is lowered (as suggested by the inflation- 
ary scenario without a cosmo logical constant) the gluino would become more 
accessible to accelerator detection. 

(iv) The largest event rates occur at large tan/? and small nig and small mo (Figs. 
7-9). Thus detectors are most sensitive to these domains. Models with very 
large tan/? (i.e. tan/3 ~50) may therefore be testable by planned detectors. 

(v) The predicted b^ s + 7 decay is large in the same region of parameter space 
(small m^, small toq) where dark matter event R is large when \i and At have 
the same sign. The current experimental rate for this decay thus eliminates 
part of the parameter space where R is large. Since the relic density constraint 
eliminates most of the A^ < part of the parameter space, the b— *> s -|- 7 
constraint signiflcantly reduces the expected event rate for ^u > 0, but does not 
effect the fx < event rates a great deal, and large event rates can still occur 
for (U < 0. 

(vi) For models possessing in addition SU(5)-type proton decay, there remain re- 
gions in parameter space, Eq. (27), satisfying both the relic density constraint 
and current proton decay bounds. These points also are within the 90% C.L. 
bounds of the current b^ s + 7 decay rate. However, since proton decay favors 
large mo and tan/3 <10, the predicted dark matter event rates are all for R 
<0.01 event/kg da. Hence models with SU(5)-type proton decay predict that 
Zi dark matter will be inaccessible to current detectors. 
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